Percolation in self-similar networks 
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We provide a simple proof that graphs in a general class of self-similar networks have zero perco- 
lation threshold. The considered self-similar networks include random scale-free graphs with given 
expected node degrees and zero clustering, scale-free graphs with finite clustering and metric struc- 
ture, growing scale-free networks, and many real networks. The proof and the derivation of the 
giant component size do not require the assumption that networks are treelike. Our results rely 
only on the observation that self-similar networks possess a hierarchy of nested subgraphs whose 
average degree grows with their depth in the hierarchy. We conjecture that this property is pivotal 
for percolation in networks. 
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Percolation is a fundamental phenomenon in nature. 
Recent developments in percolation theory [1] open new 
perspectives in many areas of statistical mechanics and 
quantum field theory [2 . In statistical mechanics of com- 
plex networks, the percolation properties of a network 
determine its robustness with respect to structural dam- 
age, and dictate how emergent phenomena depend on the 
network structure [3] . Large clusters of connected nodes 
emerge above a critical value of some network parameter, 
e.g., the average degree; below the threshold, networks 
decompose into a myriad of small components. This per- 
colation threshold can be zero, meaning that networks 
are always in the percolated phase. A classic example 
is random scale-free networks with the power-law degree 
distribution exponent 7 lying between 2 and 3 [H [5] . The 
value of the percolation threshold, the size of the giant 
component, and the specifics of the percolation transi- 
tion strongly depend on fine details of the network topol- 
ogy [3] . This dependency hinders attempts to define per- 
colation universality classes, even though some networks 
show some degree of percolation universality [B]. 

This problem is aggravated by difficulties in the an- 
alytic treatment of percolation properties for networks 
with strong clustering. A majority of the obtained ana- 
lytic results use the generating function formalism based 
on the assumption that networks are locally treelike [TJ. 
This assumption allows one to employ convenient tools 
from the theory of random branching processes. The as- 
sumed absence of loops implies, in particular, that clus- 
tering is zero in the thermodynamic limit. This zero- 
clustering approximation is valid for weakly clustered 
networks where triangles do not overlap, but it is invalid 
for networks with strong clustering and overlapping tri- 
angles observed in many real systems [5] . Noticeably, the 
exact results derived for some network models with clus- 
tering can be mapped to treelike zero-clustering graphs 
after appropriate transformations [5]. 



In this Letter, we provide a remarkably simple rigorous 
proof for the absence of a percolation threshold in a gen- 
eral class of self-similar networks. The proof does not rely 
on the treelike assumption or on generating functions. 
It does not depend on whether a network is weakly or 
strongly clustered, and it applies equally well to equilib- 
rium or non-equilibrium networks. The proof relies only 
on network self-similarity, defined as statistical invariancc 
of a hierarchy of nested subgraphs with respect to a net- 
work renormalization procedure. The percolation thresh- 
old is zero as soon as the average degree in subgraphs is a 
growing function of their depth in the hierarchy — a prop- 
erty characterizing many real networks. We also calculate 
analytically the size of the giant component, supporting 
all the results by large-scale numerical simulations. 

Let Q{{a}) be an ensemble of sparse graphs in the 
thermodynamic limit, where {a} is the set of model pa- 
rameters. In the case of classical random graphs, for ex- 
ample, set {a} is just the average degree (fc). Consider a 
transformation rule T that for each graph G E G({a>}) se- 
lects one of G's subgraphs. Denote the ensemble of these 
subgraphs by Qt{{oi})- The ensemble G({a}) is called 
self-similar with respect to the transformation rule T if 
the transformed ensemble is the same as the original one 
except for some transformation of the model parameters, 

<MM) = Q{{a T }). (1) 

In what follows we describe three general types of graphs 
to which this definition applies. The first two types are 
equilibrium random scale- free graph ensembles belonging 
to a general class of network models with hidden vari- 
ables |10j . The third one is a non-equilibrium ensemble 
of growing networks. 

Type I: The graphs in this ensemble are constructed 
by assigning to each node a hidden variable k drawn from 
the power-law probability density p{n) = (7— l)/tg k / . 
Without loss of generality, k can be selected such that 
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k > kq is the expected degree of nodes with hidden 
variable k, so that the degree distribution scales as a 
power law with exponent 7. Each pair of nodes with 
expected degrees k and k! is then connected with proba- 
bility r(«, k') = f((j,Kn'), where constant p, fixes the av- 
erage degree (k) in the constructed graphs, and function 
fix) < 1 is an arbitrary analytic function with /(0) = 0. 
This type of graphs includes as particular cases the max- 
imally random graphs with a given expected degree se- 
quence [11] , and random graphs with arbitrary structural 
correlations |TS] . In the former case, 
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Clustering vanishes in the thermodynamic limit, and 
therefore the treelike assumption holds. 

Type II: Besides having assigned expected degrees 
k, nodes in this type of graphs are also uniformly dis- 
tributed in a homogeneous and isotropic D-dimensional 
metric space [13] . Here, we consider a circle of radius 
R with a constant density of nodes 6 — N/(2%R), al- 
though all the following results can be extended to an 
arbitrary dimension. The connection probability between 
a pair of nodes with hidden variables k and k' separated 
by distance d = (n — \n — \6 — 0'\\)R on the circle (0's 
are the node angular coordinates) must be of the form 

r(re, 9; k', 6') — h (7^7^ , where function h must be inte- 

grable. These graphs have the same degree distribution 
as the type I graphs, but clustering is finite in the ther- 
modynamic limit, thanks to the triangle inequality in the 
underlying metric space [13] . Therefore the treelike as- 
sumption does not hold. 

The graph sparsity in the thermodynamic limit defines 
constant fi in the two cases as 
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where I — J h{x)dx. Since kq and <5 are dumb param- 
eters, we see that unless functions / and h contain some 
additional parameters, the described two graph ensem- 
bles have only two independent parameters: the power- 
law exponent 7 and the average degree (k). 

Consider now transformation rule T which simply re- 
moves all nodes with hidden variable k < kt from 
a given graph G in any of the two ensembles, where 
kt is some predefined threshold. This transformation 
maps the original graph G to its subgraph Gt of size 
Nt = N(kq/ kt) 7 ^ 1 ■ The hidden variables K of nodes 
remaining in Gt are distributed according to pt(k) — 
(7 — 1)k^ -1 k -7 with k > kt- That is, the power-law ex- 
ponent in Gt is the same as in G, 7t = 7- The transfor- 
mation does not affect the hidden variables of the nodes 
in subgraph Gt- Therefore the connection probability in 
Gt is exactly the same as in the original graph G, which 



means that the ensemble of transformed graphs is iden- 
tical to the ensemble of original graphs, except that the 
average degree has changed. Specifically, the transfor- 
mation of parameters {a} —> {aj-} in the self-similarity 
definition in Eq. ([I]) is 



7 -> It = 1, (k) -> (k) T = (k) 



N 



(4) 



which is the same for both type I and type II graphs [13] . 
Therefore, both ensembles belong to the same self- 
similarity universality class. 

Type III: As opposed to the first two equilibrium en- 
sembles, the graphs of this type are grown by adding 
nodes one by one. Each node i brings rrii new links, 
where rrii — mo(N /i) n and rj € [0, 1). Each link is then 
attached to a random existing node. A network is ini- 
tialized with lOmo-ZV uncounted disconnected nodes. If 
rj = 0, the generated graphs have an exponential de- 
gree distribution. If rj > 0, the degree distribution 
is P(k) = (l + m /r]) 1 / r i/[i 1 (k + mo/r)) 1+1 /% i.e., ap- 
proximately a power law with exponent 7 = 1 + 1/rj. 
The transformation rule T simply extracts from a grown 
graph its subgraph composed of the first Nt nodes. With 
this T, the graphs of this type are also self-similar, and 
the parameter transformation in definition (JlJ is 



7 -> 7T = 7, (fc) ->■ (k) T = (k) 



N 
Nt 



(5) 



As we show next, self-similarity of the considered en- 
sembles (types I, II, and III), and the proportionality 
Nt ~ N are sufficient to prove the absence of a per- 
colation threshold for equilibrium scale-free graphs with 
exponent 7 < 3 and for growing graphs with any 7. The 
key property which we will use is that the average de- 
gree of these self-similar subgraphs is a growing function 
of the subgraph depth in the nested subgraph hierarchy, 
meaning that (k) T grows as N T decreases in Eqs. (4|5 1. 
The same property characterizes many real networks as 
shown in Fig. [I] and in [BHH]. 

The proof is by contradiction. As usual 0], let the 
average degree be the order parameter for a percolation 
transition. Suppose that the considered self-similar en- 
sembles do have a non-zero percolation threshold at some 
critical value of the average degree (k) c . Consider a graph 
with the average degree below the threshold ((fc) < (k) c ) 
which has no giant component. Since its subgraphs be- 
long to the same ensemble, their percolation threshold 
is also (k) c . But since their average degree increases 
with their depth in the subgraph hierarchy, there exist 
deep enough subgraphs whose average degree is above the 
threshold ((&}t > (k) c )- We thus arrive at a contradic- 
tion since a graph which does not have a giant component 
must contain subgraphs which do have giant components. 

We next compute the size g{p) of the giant component 
in bond percolations with bond occupation probability p, 
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FIG. 1: Ratio of the subgraph average degree {k)r to the 
average degree (k) in the whole graph as a function of the in- 
verse relative subgraph size N/Nt for few real networks. The 
subgraphs are obtained by removing nodes with degrees below 
thresholds kx from the original network. To insulate against 
finite-size effects, the data is shown only for subgraphs of size 
Nt/N > 0.1. Actors, actor collaborations from the Inter- 
net Movie Database; Airports, USA airport network; English, 
web of semantic associations between words in English; In- 
ternet, topology of the Internet at the Autonomous Systems 
level; Proteins, protein interaction network of Saccharomyces 
cerevisiae; and Trust, mutual trust relationships among indi- 
viduals extracted from the Pretty Good Privacy data. 

confirming the absence of the percolation threshold in the 
considered ensembles. We first focus on the equilibrium 
networks of types 1 and If, in which case bond percolation 
is equivalent to replacing the connection probability rij 
with prij. Given a node i and a set of other nodes T, the 
probability that i is connected to at least one node in T is 
one minus the probability that i is not connected to any 



node in T, i.e., 1 — exp 



Node i be- 



EjexW 1 -Wij) 
longs to the giant component if and only if it is connected 
to the giant component of the graph without i. lfgj(p) 
denotes the probability that this i-deprived component 
contains some other node j, then 



i(p) = 1 - exp 



J29j(p) In (1 - pry) 



(6) 



Since in small-world networks a single node cannot signif- 
icantly affect the percolation properties of the rest of the 
graph, we identify cjj(p) — 9j(p), transforming Eq. ^ 
into a self-consistent equation for gi(p). We note that 
Eq. ^ does not use the treelike assumption. This equa- 
tion is thus valid for the type II graphs with strong clus- 
tering as well as for zero-clustering type I graphs. It leads 
in the thermodynamic limit to the following expression 
for the probability g(n;p) that a node with expected de- 
gree k belongs to the giant component: 
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FIG. 2: Relative size of the largest cluster {Si)/N vs. the 
analytical solution for type II networks with 7 = 2.5, (k) = 
3(7 — 1)/(t — 2), h(x) — e~ x , and average clustering coefficient 
(measured over degrees larger than k = 1) 5 = 0.5. 

for types I and II, respectively. The size of the giant 
component is then g(p) = J p{K)g{n;p)dK = 1 — (7 — 
l)i? 7 [V>(p)], where is the extended exponential inte- 
gral. In diluted networks with p <C 1, ai(p) ~ au{p), 
and the giant component size for both classes becomes 
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The value of the critical exponent (3 in g{jp) ~ p@ is thus 
P = 1/(3 — 7), agreeing with [15 . We emphasize that in 
our case, this result is obtained without using the treelike 
assumption. Therefore, quite surprisingly, this exponent 
characterizes equilibrium scale-free networks with arbi- 
trary clustering and degree correlations. 

In non-equilibrium networks of type III, we can com- 
pute an upper bound for /3. Self-similarity of these net- 
works, coupled with the observation that any node be- 
longing to the giant component of a self-similar subgraph 
of a type III graph belongs also to the giant component 
of the graph itself, leads to inequality 

M*W°{[nt\ p )- (12) 

By choosing Nt/N — p 1 ^ , we obtain g(p) > p^^/^gil). 
Therefore the exponent j3 satisfies j3 < 1/rj = 7— 1. 
We see that growth reduces significantly this exponent, 
compared to the equilibrium case with the same 7. 

We next check our analytic results against large-scale 
simulations. We generate type I and II networks using 
the connection probabilities in Eq. ^ and h(x) = e~ x , 
respectively. We do not allow k's above the natural cutoff 
k c = N 1 /^' 1 ^. For all the three graph types, for each 
graph size N ranging from 10 3 to 10 5 , and for each value 
of the bond occupation probability p, we generate 10 3 
graphs, and for each graph we perform bond percolation 
10 4 times. For each percolation we measure the size Si of 
the largest connected component in the graph using the 
fast algorithm of Newman and Ziff [16 , and calculate 
the average (Si) of the largest component size and its 
fluctuations, i.e., susceptibility \ = \j ((Si ~ (Si)) 2 ), for 
each combination of p, N, and graph type. 
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FIG. 3: Bond percolation simulations for equilibrium (types I 
and II) and growing (type III) networks . a: susceptibility \ 
as a function of bond occupation probability p and graph size 
N for the same network as in Fig. [2] b and c: position p-max 
and height ^ max of the peak of x as functions of network size 
N. The straight lines are power law fits, d: exponents 1/v* 
and 7*/^* in p max (N) ~ N" 1 ^' and X ™*(N) ~ AP*/"* 
for the type I and II graphs, e and f: Bond percolation 
simulations for non-equilibrium networks (type III) with 77 = 
1/4 (7 = 5) and mo = 2. The measured values of the scaling 
exponents are 1/v* — 0.24(3) and 7*/^* = 0.3(8). 

Since the convergence to the thermodynamic limit in 
scale-free networks is slow [T7] , it is difficult to accurately 
measure exponent f3 in simulations. Nevertheless we ob- 
serve an agreement, albeit slowly converging, between 
the analytical solution for g{p) and simulations in Fig. [2] 
In Fig. [3] we also show susceptibility x(pi N) for equilib- 
rium and growing networks. Susceptibility displays peaks 
whose positions p m ax{N) and heights x max (N) depend 
as power laws on the system size, p m ax(N) ~ N^ 1 / 1 ' 
and x max i N ) ~ N^'/ v ' . Taken together, these two 
results confirm that the giant component emerges at 
P = Pmaxi aud that the percolation threshold vanishes 
in the thermodynamic limit N — > 00 where p m ax ~> 
and x max -> 00. 

In short, self-similar networks with subgraphs of grow- 
ing average degree have no percolation threshold. The 
proof can be generalized to any processes with phase 
transitions whose critical points depend monotonously 
on the average degree. Examples include, among others, 
the absence of an epidemic threshold in epidemic spread- 
ing processes, or the absence of a paramagnetic phase in 
the Ising model on scale-free networks [3|. 

The identification of percolation universality classes 
for general random networks is a notoriously difficult 
problem — details tend to prevail. Nevertheless, the re- 



sults presented here lead us to conjecture that self-similar 
networks can be split into three general percolation uni- 
versality classes, depending only on whether the aver- 
age degree in the nested subgraph hierarchy increases, 
remains constant, or decreases with the subgraph depth, 
and independent of any other network properties, such as 
clustering, correlations, equilibrium vs. non-equilibrium 
classification, etc. 
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